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ABSTRACT 

We investigate the criteria for successful core-collapse supernova explosions by the neutrino mech- 
anism. We find that a critical-luminosity/mass-accretion-rate condition distinguishes non-exploding 
from exploding models in hydrodynamic one-dimensional (ID) and two-dimensional (2D) simulations. 
We present 95 such simulations that parametrically explore the dependence on neutrino luminosity, 
mass accretion rate, resolution, and dimensionality. While radial oscillations mediate the transition 
between ID accretion (non-exploding) and exploding simulations, the non-radial standing accretion 
shock instability characterizes 2D simulations. We find that it is useful to compare the average dwell 
time of matter in the gain region with the corresponding heating timescale, but that tracking the 
residence time distribution function of tracer particles better describes the complex flows in multi- 
dimensional simulations. Integral quantities such as the net heating rate, heating efficiency, and mass 
in the gain region decrease with time in non-exploding models, but for 2D exploding models, increase 
before, during, and after explosion. At the onset of explosion in 2D, the heating efficiency is ^2% 
to ~5% and the mass in the gain region is ~0.005 Mq to ~0.01 M0. Importantly, we find that the 
critical luminosity for explosions in 2D is ~70% of the critical luminosity required in ID. This result 
is not sensitive to resolution or whether the 2D computational domain is a quadrant or the full 180°. 
We suggest that the relaxation of the explosion condition in going from ID to 2D (and to, perhaps, 
3D) is of a general character and is not limited by the parametric nature of this study. 

instabilities — methods: numerical — shock waves — supernovae: 



Subject headings: hydrodynamics 
general 

1. INTRODUCTION 

Four decades of core-collapse simulations have in- 
creased our understanding of the core-collapse mech- 
anism, yet a complete theory of the mechanism has 
not emerged. Dete c tion of neutrinos d uring SN 1987A 
(jBionta et all Il987l : iHirata et all Il987f ) confirmed only 
the fundamentals of core-collapse supernovae. Theory 
suggests that the Fc core collapses to form a protoneu- 
tron star (PNS), which launches a shock wave. Before the 
bounce shock can explode the star, it is sapped of energy 
by nuclear dissociati on and neutrin o losses and staUs int o 
an accretion shock (|Mazureklll9"8l : lBruennlfT985lll989f ). 
Understanding the mechanism that revives the stalled 
shock into explosion has been the goal of the community 
for many decades. 

Since the pio n eering work of IWilsonI (fl98l and 
iBethe fc WilsonI (|1985f ). the favored shock revival 
mechanism has been the delayed-neutrino mechanism 
(or simply neutrino mechanism), in which neutrinos 
heat the post-shock region and restart the shock's 
outward progress after hundreds of milliseconds of delay. 
Detailed one-dimensional (ID) simulations using state- 
of-the-art equations of state (EOSs), neutrino-matter 
cross sections, and neutrino transport have shown that 
the neu trino mechanism f a ils to pr oduce explosions 
in I D (iLiebendorfer et all l2001bl[al: iRampp fc Jankal 
I2002t iBuras et al.l 120031 [Thompson et all l2003l : 



ILiebendorfer et all l2005l). except the least massive 



of the massive stars (|Kitaura et aI1l2006l : iBurrows et al.l 
l2007a( ). Recent 2-dimensional (2D) simulations, and the 
accompanying aspherical instabilities, have suggested 
that the neutrino mechanis m may ye t succeed, though 
it fa ils in ID dHcrant et"an 11994 [janka fc Miille^ 
19951 IBurrows etall 119951 iJanka fc Miilled Il996l: 
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2006bllal : iMarek fc Jankal l2007t lOtt et al. r|200^. ThusT 
the fundamental question of core-collapse theory 
remains; how is accretion reversed into explosion? 

Exposing the core-collapse mechanism will require de- 
tailed three-dimensional (3D) radiation-hydrodynamic 
simulations. However, the core-collapse problem is 
messy, with many subtle nonlinear couplings and feed- 
backs on local and global scales, and extracting the 
essence of the mechanism from even ID radiation- 
hydrodynamic simulations has proven very difficult. 
Hence, revealing the core ingredients and conditions 
of the mechanism will likely require a two-front at- 
tack. On the one hand, multi-dimensional radiation- 
hydrodynamic simulations will give fully consistent ex- 
plosions, producing observationally testable energies, 
neutron star masses, nucleosynthetic yields, and more. 
On the other hand, simplified models that nevertheless 
retain the important physics, but allow adjustment of im- 
portant parameters will help reveal what is important. In 
this paper, we pursue the latter philosophy. We present a 
systematic parameterization of the conditions and crite- 
ria for explosion by the neutrino mechanism, emphasizing 
the effect that going to 2D from ID has on the critical 
n eutrino luminosity r equire d for explosions. 

iBurrows fc Goshvl (|1993f ) suggested a simple frame- 
work for determining the conditions for successful explo- 
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sions by the neutrino mechanism. They approximated 
the stalled shock and accretion phase as a steady-state 
problem, transforming the governing partial differential 
equations into ordinary differential equations. By pa- 
rameterizing the electron-neutrino luminosity, L^^ , and 
the mass accretion rate, M, they identified a critical 
Ly^-M curve that distinguishes steady state accretion 
solutions (lower luminosities) from explosions (high lu- 
minosities). This implied that global conditions, not lo- 
cal conditions, mediate the transition from accretion to 
explosion, which in turn suggests that core-collapse ex- 
plosion is a global ins t ability. Mor e than a decade later, 
lYamasaki fc Yamadal ()2005l . l2006f ) reproduced these re- 
sults and performed a linear stability analysis of the 
steady-state solutions that showed unstable solutions 
near the critical luminosity. In this paper, we extend 
this previous work by performing a suite of hydrody- 
namic simulations in both ID and 2D. 

Assuming that the concept of a critical luminosity 
applies to 2D simulations, many have noted that crit- 
ical luminosities for 2D simulations are likely to be 
lower than in ID simulations. In their ID s teady- 
state solutions. lYamasaki fc Yamadal (|2005l . [20061 ) inves- 
tigated the effects of rotation and convection on the 
critical luminosity, but these analyses lacked a self- 
consistent treatment of what is an inherently 2D/3D phe- 
nomenon. Early 2D simulations using flux-limited diffu- 
sion noted a trend toward ex plosions aided by neutrino- 
driven convection (iBurrows ct al. 1995; J anka fc Miiller, 
Il996l) . while ID simulations failed to e xplode for the 
same neutrino luminosity. More recently, iBlondin et al.l 
(I2OOI identified a new instability that may prove cru- 
cial for the viability of the neutrino mechanism. It is 
the standing accretion shock i nstability (SASI) , which 
may be an advective-acoustic ([Foglizzo fc Tagged I2OOOI: 
lFoglizzoll2002l: iFbglizzo et al.|2007|[) or purelv acoustic cy- 
cle (jBlondin fc Mezzacappall2006[ ). Whichever explains 
the SASI, recent simulations have suggested that it may 
facilitate the neutr ino mechanism (Marek & Janka 2003; 
iBuras et al.ll2006ah . However, these simulations have yet 
to demonstrate a reliable explosion mechanism for a wide 
range of progenitor masses that yields explosion ener- 
gies co nsistent with Nature. Furthermore ^Burro ws et sl\ 
(120061). IBurrows et al l (l2007dl). [B urrows ct al. (2007cD. 
IBurrows et al.l (l2007al) . and IBurrows et alT (l2007bl) . us- 
ing VULCAN/2D (jLivnel [19931 iLivne et all 120041 ) and 
multi-group flux-limited diffusion (MGFLD) , obtain suc- 
cessful explosions only for the 8.8-M (7^ model by the 
neutrino mechanism alone. Recently, lOtt et al] (f2008l ) 
compared MGFLD and multi-angle transport in VUL- 
CAN/2D. While they note some interesting differences 
between the transport schemes, explosions remain elu- 
sive. Moreover, it has not been demonstrated that the 
concept of a critical L^^-M condition pertains to 2D sim- 
ulations, and if it does, how do the critical luminosities 
of ID and 2D simulations compare? 

To address these questions, we conduct ID and 2D 
simulations for various values of L^^. In the past, there 
have been a few investigations on the systematic ef- 
fects of neutrino lumin osity on the explosio n mechanism 
(iJanka fc Miillerl l l996l). neutron star kicks (jScheck et al.l 
120061 ). and the SASI (johnishi et al.ll2006( ). but none has 
thoroughly investigated both L^^ and M in ID and 2D 



simulations to address the viability of a critical luminos- 
ity condition for explosion. In addition to and M, 
we compare simulations with differe nt spatial re solutions. 
Using the code, BETHE-hydro (jMurphv fc Burrowi 
[2008I) . we simulate the core-collapse, bounce, and post- 
bounce phases in time-dependent ID and 2D simu- 
lations. These simulations have no inner boundary 
and include the PNS core. A finite-temperature EOS 
that accounts for nucleons, nuclei, photons, electrons, 
posit rons, and all the appropriate phase transitions is 
used (|Shen et al.ll998l). Empl o ying 11.2- and 15-Mp pro- 
genitors (jWooslev et al.ll2002l : [Wooslev fc Weaveilll995f l 
as initial conditions, a wide range of M is sampled (from 
--0.08 Mq/s to ~0.3 Mg/s). Finally, we use stan- 
dard approximations for neutrino heating and cooling 
that enable a straightforward par ameterization of L^^ 
(iBethe fc Wilsonl[l985t lJankal[200ll) . 

The basic equations and the numerical techniques are 
presented in ^ In ^ we discuss the progenitor mod- 
els and describe the suite of simulations performed. In 
21 we analyze the evolution of the shock for ID and 2D 
simulations and discuss the role shock oscillations play 
in the explosion. In ^ we investigate the effect that the 
grid, specifically spatial resolution and the angular ex- 
tent of the domain, has on explosions. In ^ we revisit 
the condition for explosion as expressed by the heating 
and advection timescales, and in SJ7| describe other con- 
ditions at explosion. In fj8l we quantify the differences in 
the critical luminosity condition between ID and 2D sim- 
ulations. Finally, in SJSJ we summarize our conclusions. 

2. EQUATIONS AND NUMERICAL TECHNIQUES 

The basic equations of hydrodynamics are the conser- 
vation of mass, momentum, and energy: 

dp 

'dt 



dv ^ ^ 

p— = -pV$- VP, 
dt 



and 



p— = -PV • V- 
dt 



n-c. 



(1) 
(2) 

(3) 



p is the mass density, v is the fluid velocity, $ is the 
gravitational potential, P is the isotropic pressure, e is 
the specific internal energy, and d/dt = d/dt + v-V is the 
Lagrangian time derivative. In this work, the neutrino 
heating, 7i, and cooling, C, terms in eq. ([3]) are assumed 
to be 
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Note that these approximations for heating and c ool- 
ing by neutrinos (|Bethe fc Wilsonlll985l : lJankall2001h de- 
pend upon local quantities and predefined parameters. 
They are p, temperature (T), the distance from the cen- 
ter (r) , the electron-neutrino temperature (T^^ ) , and the 
electron- neutrino luminosity (L^^), which is in units of 
Iq52 g]-g g-i gy using eqs. ^ and (O, we gain consider- 
able time savings by approximating the effects of detailed 
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neutrino transport. For all simulations, we set T^^ ~ 4 
MeV. In eq. (0]), it has been assumed that L^,^ = Ly^ 
and that the mass fractions of protons and neutrons sum 
to one. Therefore, the sum of the electron- and anti- 
electron-neutrino luminosities is ii/^p^ = 2Lu^. Closure 
for eqs. (IT][3]) is obtained with an EOS appropriate for 
matter in or near nuclear statistical equilibrium (NSE) 
(jShen et al.lll998[ ). and the effects of photons, electrons, 
and positrons are included. As such, the EOS has the 
following dependencies: 

F = P(p,£,re), (6) 

where is the electron fraction. Therefore, we also solve 
the equation: 



where Te is the net rate of Ye change. 

Using BETHE-hydro (|Murphv fc Burrows! [2001 . we 
solve eqs. ([T][3l) in one and two dimensions by the Ar- 
bitrary Lagrangian-Eulerian (ALE) method. To advance 
the discrete equations of hydrodynamics by one timestep, 
ALE methods generally use two operations, a Lagrangian 
hydrodynamic step followed by a remap. The structure 
of BETHE-hydro's hydrodynamic solver is designed for 
arbitrary-unstructured grids, and the remapping com- 
ponent offers control of the time evolution of the grid. 
Taken together, these features enable the use of time- 
dependent arbitrary grids to avoid some unwanted fea- 
tures of traditional grids. 

For the calculations presented in this paper, we use this 
flexibility to avoid the singularity of spherical grids in two 
dimensions. While spherical grids are generally useful for 
core-collapse simulations, the convergence of grid lines 
near the center place extreme constraints on the timestep 
via the Courant-Friedrichs-Levy condition. A common 
remedy is to simulate the inner '--^lO km in ID or to use 
an inner boundary condition. Another appr oach, which 
has been used in VULCAN/2D simulations (|Livndll993l : 
iLivne et"ani2004f ). is to avoid the singularity with a grid 
that is pseudo-Cartesian near the center and smoothly 
transitions to a spherical grid at l arger radii. We use 
a sim ilar grid, the butterfly mesh (jMurphv fc Burrows! 
|2008|) . for the simulations in this paper. 

For ID simulations, 700 radial zones are distributed 
from the center to 4000 km. The innermost 100 zones 
have a resolution of 0.34km, and the remaining 600 zones 
are spaced logarithmically from 34 km to 4000 km. Af- 
ter each Lagrangian hydrodynamic solve, the flow is 
remapped back to the original grid, which in effect pro- 
duces an Eulerian calculation. For 2D simulations, we 
use a butterfly mesh interior to 34 km and a spherical 
grid exterior to this radius. For all 180° simulations, the 
innermost pseudo-rectangular region is 50 by 100 zones, 
and the region that transitions from Cartesian to spher- 
ical geometry has 50 radial zones and 200 angular zones. 
The outermost spherical region has 200 angular zones 
and 300 or 150 radial zones which are logarithmically 
spaced between 34 km and either 4000 km or 1000 km. 
The butterfly portion has an effective radial resolution 
of 0.34 km with the shortest cell edge being 0.28 km. 
The grids for the 90° simulations are similar, but with 
half the number of zones. In 2D simulations, we do not 
remap every timestep. Since the timestep is limited by 
the large sound speeds in the PNS core, we can afford 



to perform several Lagrangian hydrodynamic solves be- 
fore remapping back to the original grid. By remapping 
seldomly, we gain considerable savings in computational 
time. 

We investigate the systematics of the neutrino mecha- 
nism by parameterizing L^^ , M, resolution, and the di- 
mensions of the simulation. Since this requires a large 
set of simulations, we make several approximations to 
expedite the calculations. For one, gravity is calculated 
via g = —GMint/r'^: where Mint is the mass interior to 
the radius r. During collapse, rather than solving the 
rate equations for Yg due to ne utrino transp o rt, we use 
a Ye vs. p parameterization. iLiebendorfeH l|2005f ) ob- 
served that ID simulations including neutrino transport 
produce Ye values during collapse which arc essentially a 
function of density alone. This allows for a parameteriza- 
tion of Ye as a fun ction of p. To change Ye , we use results 
of ID SESAME (jBurrows et all l2000l : [Thompson et all 
l2003l ) simulations to defin e the function Yp(p \ and we 
employ the prescription of ILiebendorfeH ([2005') to calcu- 
late local values of Fg. As a result, the most important 
effects of electron capture during collapse are included 
without the need for expensive, detailed neutrino trans- 
port. 

3. PROGENITOR AND SIMULATION MODELS 

To sample a range of M, we initiate the simulations 
with two core-collapse progenitor models that have dif- 
ferent density profiles: a 15 M0 proge nitor with a shallow 
density profile exterior to the Fe core ijWooslev fc Weaveil 
Il995l ). and a 11.2 Mf^ progenitor with a steeper exte- 
rior density profile (jWooslev et al.|[2003 ). In Fig. [U the 
mass accretion rate, M, vs. time at 250 km (just above 
the stalled shock) is shown for ID non-exploding models. 
The solid and dashed lines show the time-dependent ac- 
cretion rate for the 11.2 and 15.0 models. While the 
outer portions of the Fe core accrete through the shock, 
M is as high as 10 Mq/s and decreases to 2 Mq/s within 
50 ms of bounce for both masses. After the Fe core is 
fully accreted, the accretion rates for the two models de- 
viate. It takes 50 ms to fully accrete the Si-burning shell 
using the II.2-M0 model and 100 ms using the 15-M0 
model. Afterward, M declines slowly to 0.2 Mq/s (0.08 
Mq/s for the 15 M© model). Together, these two models 
slowly sweep through a range of accretion rates from 0.3 
Mq/s to 0.08 Mq/s, enabling a study of the M depen- 
dence. 

The results presented in the following sections are de- 
rived from 95 simulations that as a group, represent a 
parameterization of L^^, M, resolution, and dimension- 
ality. Table [1] lists these simulations into 14 sequences. 
Each sequence of simulations is distinguished by the pro- 
genitor model (column 2), dimensionality (column 3), 
number of radial zones (column 4), radius of the outer 
boundary (column 5), and a range of electron- neutrino 
luminosities (column 6). The sequence names encode the 
information contained in columns 2 through 5. The first 
four characters indicate the progenitor model. Next are 
the dimensionality labels: "ID" for ID simulations, "2D" 
for 2D simulations that incorporate the full of 6 (180°), 
and "Q" for 2D simulations that simulate the quadrant 
that lies between the pole and equator (90°). For the 
2D simulations, the final character represents the resolu- 
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tion: "1" for the fowcst rcsolutfon of 250 effective radial 
zones within 4000 km, "2" for a finer resolution of 400 
effective radial zones in 4000 km, and "3" for the highest 
resolution of 400 effective radial zones within 1000 km. 
On the whole. Table [T] represents a parameter study for 
the conditions near explosion in L^^^ M, dimensionality, 
and resolution. 

4. RADIAL SHOCK OSCILLATIONS IN ID AND THE SASI 

IN 2D 

When iBurrows fc Goshvl (|1993f ) reported a critical lu- 
minosity vs. accretion rate condition for successful ex- 
plosions, they did so based upon steady-state solutions. 
It was left to subsequent calculations to show that the 
concept of a critical luminosity condition was appropri- 
ate for time-dependent simulations. With plots of the 
evolution of the shock radius, we show that a critical lu- 
minosity condition distinguishes non-exploding and ex- 
ploding models in time-dependent ID and 2D simula- 
tions, and that the critical luminosity for 2D simulations 
is lower than the critical luminosity for ID simulations. 
In addition, we report global oscillations near the critical 
luminosity in both ID and 2D simulations, but we show 
that the oscillations are of a quite different character. 

The time evolution for the shock radius, i?shock, of ID 
simulations is presented in Fig. [2l In the top panel, 
we display the radii for many models of the 15.0-lD se- 
quence, and each is labeled by L^^ . The range of lu- 
minosities shown, L^^ = 2.2-2.9 (in units of 10^^ erg 
s~^), highlight values near the critical luminosity for 
this model. The fact that several luminosities lead to 
explosion at different times is an early indication that 
the critical luminosity depends upon the mass accretion 
rate. In f|8l we elaborate on the accretion rate depen- 
dence. Similarly, the bottom panel shows i?shock for the 
range of luminosities surrounding the critical luminosity, 
L^^ = 1.1-1.7, of the 11.2-lD sequence. 

Prior to the accretion of the Si/0 interface, the shock 
stalls for both sequences at radii that are slightly de- 
pendent upon L^^ . The accretion of this interface and 
the abrupt drop in M trigger an immediate explosion 
for the highest luminosity, = 1.7 (2.9) for the 11.2- 
ID (15.0-lD) sequence and osciUations in i?shock for the 
other luminosities considered. These radial oscillations 
either decay or maintain large amplitudes until the model 
explodes. As an example, the oscillation periods of the 
15.0-lD sequence range from ^--^90 ms (L^^ = 2.2) to ^170 
ms {Li,^ = 2.8), and the timescales for decay range from 
450 ms\L^^ = 2.2) to 1.0 s (L^^ = 2.5). These radial os- 
cillations and their large amplitudes near explosion sup- 
port, but do not prove, the hypothesis that a global insta- 
bility is responsible for the transition between accretion 
and explosion. 

Near the critical luminosity, others have reported large 
radial oscillations in ID simulations. lOhnishi et al.l 
(|2006f ) executed time-depe ndent ID simulations of the 
steady state solutions in lYamasaki fc Yamadal ()2005f ) 
and observed similar pulsations, with amplitudes of or- 
der 10%. For some luminosities, we report similar ampli- 
tudes (Fig. [21), but just as many models have large am- 
plitudes of order unity. Furthermore, we find, as do they, 
that the amplitude s, periods, and L^^ arc correlated. 
iBuras et al.l (|2006bf ) simulated the collapse, bounce, and 
postbounce phases for the 15 Mq model using a Boltz- 



mann transport algorithm in ID and a "ray-by-ray" 
derivative of the same transport algorithm for 2D simula- 
tions. Generically, they did not obtain explosions. How- 
ever, for a few runs, they omitted the velocity-dependent 
terms in their transport formulation, and this resulted 
in explosions. Explosions occurred soon after the shock 
stalled in their 2D simulations, aborting any obvious os- 
cillations in shock radius. On the other hand, like our 
simulations, accretion of the Si/0 interface initiates pul- 
sations with amplitudes ^^25% in their ID simulations 
oscillations. . 

In reporting the temporal morphology of shocks in 
2D simulations, we decompose the shock position, 
Rshocki(^ , t) , into spherical harmonics: 

Rshock{e,t) = V4^^af(t)r,o , (8) 

e 

where 

aeit) = [ RsbockiO,t)Yeodn, (9) 
V47r J 

are the time-dependent coefficients, and the normaliza- 
tion for ag in eq. ([9]) ensures that oq is the average shock 
radius. As is usual, Yio are the spherical harmonics, az- 
imuthal symmetry dictates to = 0, and the spherical har- 
monics arc normalized to satisfy 

J Y^lYeodn = 1 . (10) 

These temporal coefficients for the 15.0-2D3 sequence 
are plotted in Fig. [31 The average shock radius, ao or 
(^shock), is shown in the top panel, and a dot indicates 
the time of explosion, texp for the exploding models. Ac- 
companying are error bars indicating an estimated un- 
certainty of ioxp (see fJTlfor the definition of tcxp)- The 
1=1 and £ = 2 components normalized by oq are shown 
in the middle and bottom panels. First of all, note that 
the range of luminosities bracketing the critical luminos- 
ity for the 2D simulations in Fig. [31 are lower than the 
range for ID simulations in Fig. [2l Further compari- 
son indicates that, unlike ID simulations, 2D simulations 
show very little oscillation in (i?shock) prior to explosion. 
In fact, oscillations in (i?shock) are miniscule for non- 
exploding models and only reach ~10% around tcxp for 
exploding simulations. 

Instead, SASI originated oscillations in ^ = 1 and £ = 2 
components dominate the non-exploding and exploding 
models. The amplitudes of the dipole component, ai/ao, 
for Li,^ ~ 1.5 remain low at ^2.5%, from bounce until 
500 ms. Around 500 ms, the amplitude grows to ~10%. 
On the other hand, models with L^^ = 1.6-1.9 show a 
steady rise to ^10% from bounce. Exploding models, for 
which Lj,_. = 1.8, 1.9, and 2.0, show an abrupt increase 
in amplitude to ~20%. For non-exploding models, ai/ao 
saturates at ~10%, and for L^_. = 1.5 and 1.6 it decreases 
to ~5%. At early times, simulations with L^^ = 1.5 
and L[y^ = 2.0 show the lowest and highest SASI am- 
plitudes, respectively, and are easily distinguished from 
the other simulations. While this would suggest a cor- 
relation in SASI amplitude with L^,^ , the other models 
do not show a similar monotonic trend. Rather their 
amplitudes occupy the region between L^,^ = 1.5 and 
L^^ = 2.0 but show no discernible trend with luminos- 
ity. The quadrupolc term, a-ij a^, is generally below ^5% 
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during non-exploding phases. For times just before ex- 
pfosion and afterward, 02/00 can reach ^^10%. 

Two important results are apparent in Figs. [5] and 
[3l A critical luminosity separates exploding and non- 
exploding models in time-dependent ID and 2D simula- 
tions, and the critical luminosity for 2D simulations is 
~70% of the critical luminosity for ID models. More ev- 
idence will be presented in fj8]that the critical luminosity 
depends upon M . For ID simulations, we observe pulsa- 
tions near the critical luminosity, but in 2D simulations, 
radial oscillations are all but absent. Instead, non-radial 
SASI oscillations dominate. 

5. EFFECTS OF THE GRID 

The conclusions of the previous section might depend 
upon the grid structure used in 2D simulations. To ad- 
dress this concern, we investigate the dependence of these 
results on resolution and the range of polar angle, 9, 
included in the simulations. Encoded in the sequence 
names (Table [T]) arc the types of grids used for 2D sim- 
ulations. For calculations that encompass the full range 
of 9 from pole to pole (180°) 9 is divided into 200 zones, 
and for simulations that include only the region between 
a pole and the equator (90°) 9 is divided into 100 zones. 
Sequences that use the 180° grid have '2D' in their name, 
and for sequences that use the 90° grid, 'Q' is in its 
stead. The numbers after '2D' or 'Q' indicate the resolu- 
tion. The lowest resolution, denoted by '1', extends the 
spherical grid from the outer edge of the butterfly mesh, 
34 km, to 4000 km with 150 radial zones. The middle 
resolution. '2', divides this same space into 300 radial 
zones, and the highest resolution, '3', divides the radii 
between 34 km and 1000 km into 300 zones. For all res- 
olutions, the space between the butterfly mesh and the 
outer boundaries are positioned logarithmically. These 
simulations do show nontrivial differences in the SASI 
and post-shock flow. However, the conclusions of the 
previous section are insensitive to changes in resolution 
and the range of 9 simulated. 

For these three resolutions. Fig. [4] shows (i?shock), 
oi/oo, and a2/ao vs. time for a single luminosity, L^^ = 
1.9, and the IS-M© progenitor model. Their general evo- 
lution is similar to what was discussed in 21 ^-iid shown 
in Fig. [3] However, there are some nontrivial differences. 

After the Si/0 interface is accreted, (i?shock) of the 
lowest resolution model (15.0-2D1) secularly creeps out 
to ~220 km over 100 ms. Then, outward progression 
stalls again, and the model doesn't explode until 581 ms 
after bounce. Both the 15.0-2D2 and 15.0-2D3 models 
momentarily stall at ~200 km, a smaller radius. 150 ms 
later, the average shock radius of the 15.0-2D2 model 
reaches ~220 km on its way to explosion. At an even 
later time, 550 ms, the shock radius of the 15.0-2D3 
model reaches ~220 km. By this measure, the trend with 
resolution is monotonic. However, this does not translate 
to monotonicity in explosion time with resolution. From 
the lowest resolution (15.0-2D1) to the highest resolution 
(15.0-2D3), the explosion times are 581, 359, and 649 ms. 
Yet, considering another luminosity, L,^^ = 2.0 (see Table 
[2) does present a monotonic trend toward later explosion 
times with higher resolution: from the lowest to the high- 
est resolution, ioxp = 325, 364. and 395 ms. From this, 
one might conclude that higher-resolution grids forestall 
explosion, yet only the highest resolution shows an ex- 



plosion for L^^ = 1.8. 

The £ = 1 and £ = 2 (the middle and bottom pan- 
els of Fig. m spherical harmonic coefficients, like the 
monopole term, show some differences with resolution, 
but present very few clear trends. Any trends that do 
appear are subtle. In the earliest phase, before accretion 
of the Si/0 interface, the highest resolution simulations 
exhibit the greatest amplitude, though all three resolu- 
tions have amplitudes in the i = 1 coefficient that are 
less than 5%. At the same time, the t ~ 2 component 
is comparable for all three resolutions. After accretion 
of the Si/0 interface, the nonradial components begin to 
rise for all three resolutions. From ~200 to ~325 ms, 
the 15.0-2D1 model shows the greatest increase in both 
nonradial components, while for 15.0-2D2 and 15.0-2D3 
they are comparable in growth. At ^325 ms, the 15.0- 
2D2 model shows a significant increase in the nonradial 
components that coincides with explosion. The nonradial 
components for the other two simulations continue a sec- 
ular increase in amplitude until their explosions at later 
times. All the while, the lowest resolution run consis- 
tently maintains a slightly larger amplitude in the £ = 1 
component, but the £ = 2 components are quite sim- 
ilar. Generally, before accretion of the Si/0 interface, 
oscillation amplitudes for oi are mildly correlated with 
resolution, and afterward it is mildly anti-correlated with 
resolution. 

Figure [5] compares 180° (15.0-2D3, solid) and 90° 
(15.0-Q3, dashed) simulations. (-Rshock) and 02/00 vs. 
time are plotted for two electron-neutrino luminosities 
L,j^ = 1.8 (purple curves) and L^^ = 1.9 (green curves). 
Once again dots in the top panel mark icxp! interestingly, 
the 90° simulations consistently explode before the 180° 
simulations. Other than for times near ~300 to ^400 
ms, all models show similar evolution in the amplitude of 
0-2 1 0.0- The one exception appears around ^300 to ~400 
ms, when the 90°, L^^ =1.9 model starts exploding and 
shows the largest amplitude. Even though the 90° sim- 
ulations explode earlier than the 180° simulations, the 
critical luminosities for the 180° and 90° simulations are 
very similar (see Table [2] and . 

Using a " ray-by- ray" neutrino transport algorithm, 
iBuras et al.l ()2006a[ ) performed simulations of the 11.2- 
Mq progenitor that teetered on the verge of explosion 
for a 90° calculation, but exploded using 180°. They 
reported no explosion with a 90° simulation but, with 
all else being equal, obtained an explosion with a 180° 
simulation. Differences between our simulations and 
theirs include neutrino transport, EOS, and hydrody- 
namic solver. While any of these could account for the 
apparent discrepancy, we emphasiz e the difference in the 
90° grids. The 90° simulations of iBuras et all (|2006af) 
included a wedge that extends plus and minus 45° of 
the equator and suppressed i = 1 and I ~ 2 modes, the 
two prominent modes of the SASI. Our 90° simulations 
range from the pole to t he equator and inhi bit only the 
£ = 1 mode. Regardless. IBuras et al.l ()2006af ) noted that 
their 90° simulation was close to explosion and that ex- 
tending the grid to 180° simply tipped the scales toward 
explosion. U sing the sam e hydro dynamics and neutrino 
transport as iBuras et al.l (j2006a[ ). but a 9 0° gr id that 
goes from the pole to the equator, iMarekl (j2007l ) found 
that 90° simulations explode insignificantly earlier than 
180° simulations. In either case, these results are con- 
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sistciit with our finding that the critical luminosity for 
simulations using 90° and 180° are very similar. 

Resolution and the angular size of the grid non-trivially 
affect the shock position and morphology in 2D simula- 
tions. For one, the lowest resolution simulations consis- 
tently maintain a slightly larger amplitude in the £ = I 
component. At early times, the lowest resolution simu- 
lations reach ~220 km before higher resolution simula- 
tions. This does not lead to a monotonic correlation in 
explosion times, though. In changing the angular size 
of the grid, the 90° simulations consistently explode at 
earlier times compared to the 180° simulations. Despite 
differences that result from grid changes, the conclusions 
of the previous section do not change. Specifically, the 
critical luminosity for the 90° and 180° simulations are 
very similar. 

6. HEATING AND ADVECTION TIMESCALES 

Measures that are related to the critical L^^ and M 
condition (jBurrows fc Goshvl fl99l are the timescale 
for matter to traverse the gain region, the advection 
timescale, Tadv, and the heating tirnescale of the mat- 
ter in the gain region, Tq. iThompsonl (|2000[ ) argued that 
ID core-collapse simulations do not explode because Tg 
is longer than Tadv In other words, core-collapse simula- 
tions fail because Tadv/^q < 1- He defined, ^ aT^/Q, 
and the advection t ime as the flow time across a scale 
height. iJankal ()2001l ) explored similar integral concepts, 
and proposed conditions for explosion that are akin to 

Tadv/Tg > 1 , (11) 

the opposite of t h e IThompsonl l|2000f ) failure condition. 
iThompson et al.l (j2005l ) applied this condition to ID 
core-collapse simulations and demonstrated that this ra- 
tio does indeed help to distinguish non-exploding mod- 
els from exploding models. They defined the advec- 
tion timescale as H/vr, where H is the pressure scale 
height, and the heating timescale as {P/p)/q, where q is 
the local net heating rate. Since, others have explored 
Tadv an d Tq as a condition for explosion in 2 D simu- 
lations (|Buras et all l2006bl : IScheck et all I2008D . While 
Tudv/Tq > 1 is to some extent useful in identifying explo- 
sions, our results suggest that it is only useful as a rough 
diagnostic. 

We define the heating time as a characteristic time for 
neutrinos to change the thermal energy of the matter in 
the gain region, where heating dominates cooling. 



For spherically symmetric flow, the advection timescale 



dr 

Tadv = / 



Q 

where Eq is the zero point energy for the EOS and 



/ {n- C)dm 



(12) 



(13) 



is the net h eating rate and is in tegra ted over the gain re- 
gion. While lScheck et all (|200l and lBuras et all (|2006al) 
incl ude l/2i;^-H$ in th e int egrand of eg . (IT^ , we su ggest , 



like lThompsCT] (|2000[) and lThompson et all (|2005D . that 
comparing the heating rate to the thermal energy eq. 
p2p is a more natural definition of a heating timescale. 
To be clear though, these differences in definition alter 
the absolute scale, but not the general trends. 



(14) 



where the the integral is integrated from the shock po- 
sition to the gain radius. For 2D simulations, many in- 
carnations of eq. have appeared in the literature. 

IScheck et al.l (I2008D use a solid-angle averaged version of 
eq. mD for t.h,,. iBuras et al.l lj2006af) use eq. ^ for the 
advection timescale in ID, but for 2D simulations they 
redefine this timescale as the time for a mass shell to tra- 
verse the gain region. We've explored three measures of 
the advection timescale. The first replaces the denom- 
inator in eq. (|14p with the solid-angle averaged veloc- 
ity, {vr) ~ {J Vrdfl) /Att. For the second, if one assumes 
steady-state fiow, then = Mgi^i^/ M is a useful mea- 
sure of the accretion timescale, where Mgain is the mass 
in the gain region and M is the mass fiux in or out of the 
gain region. For simplicity, we use AI just exterior to the 
shock. The advantage of this definition is that it is rel- 
atively straightforward to measure Mgain and M for 2D 
simulations. Prior to explosion, these two definitions give 
similar results, but as the model explodes, the assump- 
tion of steady-state accretion is violated and r,„ begins 
to underestimate Tadv Neither of these definitions, how- 
ever, fully capture the complex multi-dimensional nature 
of the post-shock flow. 

We define a new timescale that better captures the 
multi-dimensional nature of the post-shock flow. Post- 
processing the simulation data, we integrate the paths of 
tracer particles for 150 ms, record their trajectories, and 
define a residence time, r^, which is the duration of each 
particle's time spent in the gain region. 50,000 particles 
are initiated at either 400 km, a radius well outside the 
shock, or 150 km, the middle of the gain region. At 
either radius, the particles are randomly distributed in 
/X = cos 9. For most situations, the particles traverse the 
gain region within the integration time, 150 ms, and have 
accreted onto the PNS. Every 10 ms, a new generation of 
50,000 particles is generated, providing a time-dependent 
distribution of Tr- 

The ratio Tadv/''"? vs. time is shown in Fig. [6] for the 
15.0-lD sequence. The top panel shows the models that 
do not explode by 1.3 s after bounce, and the bottom 
panel shows those that do explode. Despite dramatic os- 
cillations, Tadv/iq rarely reaches a value of one for the 
non-exploding models (top panel). In all non-exploding 
cases, except L^^ = 2.5, the oscillations in this ratio de- 
cay with a timescale commensurate with the decay times 
of the shock radius oscillations (Fig. |2|). For exploding 
models (bottom panel), Tadv/Tq executes large excursions 
from ~0.1 to ^100. Note that Tadv/Tq makes several 
excursions to values above order unity and back again. 
Smaller luminosities produce more excursions, with one 
at Ljx^ = 2.9 and four for L^_. ~ 2.6. 

The ratio of these timescales for 2D simulations (15.0- 
2D3 sequence) are shown in Fig. [71 The models that 
don't explode, L,^^ ~ 1.5, 1.6, and 1.7, show relatively 
constant ratios for all time. On the other hand, the ex- 
ploding models, = 1.8, 1.9, and 2.0, start at a low 
ratio (~0.1) and secularly increase to ratios of order ^10. 
Unlike the ID models, once this ratio reaches large val- 
ues (in an average sense), it remains there for the rest of 
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the simulation. 

Generally, the condition Ti^dv/Tq > 1 is a useful diag- 
nostic of explosion, but given the ambiguities in defining 
Tadv and Tq, its accuracy is limited. For example. ID sim- 
ulations that explode make dramatic excursions above 
and below this line before explosion. It might be more 
sensible to state that this condition should be satisfied 
for a finite time before explosion occurs. However, one 
shouldn't expect rigor when defining a heating timescale. 
Note that this condition assumes a quasi-steady state, 
and the large amplitude oscillations in ID complicate 
this assumption. Finally, T^^dv/Tq, in 2D exploding simu- 
lations, increases monotonically for ^100 ms, while this 
ratio remains low for 2D non-exploding simulations. This 
distinction implies that 2D exploding models are well on 
their way to explosion before r^dv/'Tq > 1 is satisfied. 
Hence, Ti^dv/Tq is a useful diagnostic for explosion, but 
may not be a rigorous condition for explosion. 

In Fig. [HI we plot for all generations of tracer particles 
the mean residence times, (r^) (top panel), and the frac- 
tion (fr) of particles with > 40 ms (bottom panel) as 
a function of time for the L^^ = 1.9 model of the 15.0- 
2D3 sequence. The particles with i?start = 400 km (150 
km) are shown with solid (dashed) lines. Both (r^) and 
fr show that the particles with i?start = 400 km and 150 
km track two fundamentally different flows. Generally, 
the particles with i?start ~ 150 km have longer average 
residence times and more particles with > 40 ms than 
the generation of particles with i?start = 400 km, which 
suggests that most of the mass that accretes through the 
shock quickly finds its way to the PNS through narrow, 
low entropy plumes. Despite the quick accretion of new 
matter, the long of particles with i?start = 150 km 
indicate that a significant amount of mass in the gain 
region does linger, resulting in more heating and higher 
entropies. 

To verify this, we compute the total amount of mass 
with high entropies and compare this to the total mass 
in the gain region. To define high entropies, we choose 
a minimum entropy of 14 ke/baryon, which corresponds 
to the maximum entropy for the ID simulation of the 
same neutrino luminosity (see Fig. I13p and gives a rough 
boundary between the regions with short (low entropy) 
and long (high entropy) residence times. The ratio of 
mass with high entropy to the total mass in the gain 
region grows from 0% to 70% between 0.2 s and 0.5 s 
after bounce. For the rest of the simulation the ratio 
remains at this level. Therefore, even though most of the 
accreted mass advects quickly through the gain region 
via the down-flowing plumes, a signiflcant fraction of the 
mass in the gain region is better characterized with long 
residence times and high entropies. 

In addition, (r^) and fr show four distinct phases, with 
each successive phase having longer (rr) and higher fr- 
Distributions of that represent these four phases are 
shown in Fig. [9l and spatial colormaps of entropy are 
shown in Figs. [TOl HH and [T^] that correspond to the 
beginning (Fig. (TUl), middle (Fig. [TT]), and end (Fig. [T^ 
of each generation's integration time of 150 ms. The top 
panel of Fig. [9]shows the distribution for four generations 
with i?start = 400 km and at t = 0.130, 0.370, 0.530, 
and 0.6 s after bounce. The bottom panel shows similar 
residence-time distributions for generations with i?start = 
150 km. 



The entropy maps in the top-left panels of Figs. [TOl 
nn and [T2] show that there is very little convection or 
SASI for the first phase {t = 0.130 s), and the distribu- 
tions of Tf (Fig. [9]) show that all particles traverse the 
gain region quickly, within ^--^9 ms. At ^^150 ms, when 
the Si/0 interface is accreted, both (r^) and fr show a 
dramatic increase. The mean residence time increases 
from ^10 ms to ~40 ms for i?start = 400 and 150 km, 
and fr increases from 0% to 35% for i?start = 400 km 
and ^--^50% for i?start = 150 km. The second phase, from 
~200 ms to ^450 ms, corresponds to the growth and 
nonlinear saturation of the SASI (top-right panels of Fig. 
[lUlini). During this phase (t = 0.370 s in Fig. O, (r^) 
for i?start = 400 km (150 km) is -22 (-32), and fr is 
-10% (-30%). The third phase, from -0.45 to -0.6 s, 
shows a rise in both (r^) and fr- The distributions in 
Fig. [S]that represent this phase, t = 0.530 s, have mean 
residence times of —50 ms (—30 ms) for i?start = 150 km 
(400 km) and fr of 50% (20%). The fourth and final 
phase is marked by a significant increase of (tics) and 
/res for -Rstart = 150 km particles and continued rise for 
^start = 400 km. In Fig. [9l the example distribution is 
shown at i = 0.6 s. For i?start = 400 km, (rres) ~ 40 ms 
and /res — 25%, while (rres) — 90 ms and /res ^ 70%. 
This phase is also characterized by a large bump of par- 
ticles with very long residence times in Fig. [51 which is 
much larger for i?start = 150 km than for i?start = 400 
km. The appearance of the bump at large Tr is one of the 
most pronounced indicators that the star has reached an 
explosive situation. 

In Fig. [ini we compare ID entropy profiles of the 15.0- 
ID, L^^ = 1.9 simulation (solid red lines) with the en- 
tropy distributions shown in Fig. 1101 For each zone of 
the 2D simulation, points mark the entropy and radius 
from the center. Generally, both ID and 2D profiles show 
the shock as an abrupt jump in entropy of the accreting 
matter, the gain region as indicated by the highest en- 
tropies, and the cooling region where entropy is signifi- 
cantly reduced as matter continues to accrete onto the 
PNS. Where ID profiles show a negative entropy gradi- 
ent, 2D simulations are convectivc. The most striking 
differences, though, are that when the SASI dominates 
the post-shock flow, t = 0.370, 0.530, and 0.600 s, the 
shock is asymmetric, and the gain region has a large 
range of entropies, in which high entropies correspond 
to regions of long residence times, and low entropies cor- 
respond to plumes that funnel matter toward the PNS. 

Other than passage through shocks, the entropy is de- 
termined by neutrino heating and cooling. Therefore, 
entropy is a good measure of the integrated history of 
net heating. Considering mass- and volume-weighted 
averages, the entropy for 2D simulations is a couple 
fcs/baryon (less than 10%) higher than ID simulations. 
However, Fig. [131 shows that 2D simulations have re- 
gions with —30% higher entropies than ID simulations 
for the same L^^ and M. In other words, the long resi- 
dence times that accompany the SASI in 2D simulations 
result in regions with higher integrated net heating and 
offer an explanation for the lower critical luminosity of 
2D simulations. 

7. CONDITIONS AT EXPLOSION 
As relevant indicators of the core-collapse mechanism, 
we investigate the net heating rate (Q), heating efficiency 
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(Q/Li^^pJ, and the mass in the gain region (A/gain)- In 
general, these quantities show distinctly different evolu- 
tions for ID and 2D simulations near their respective 
critical luminosities. Plots of the ID simulations reflect 
the very large radial shock oscillations, while plots of 
the 2D simulations show relatively smooth evolution and 
bifurcations distin guishing expl o sive m odels from non- 
explosive models. IScheck et al.l (|2008[ ) compared these 
quantities for ID and 2D simulations with the same neu- 
trino luminosity and found similar bifurcations in these 
integral measures. We show that this bifurcation is 
generic when we include ID and 2D non-exploding and 
2D exploding models. In Table [21 we summarize the con- 
ditions at explosion, but first we take care to define iexp, 
the time of explosion. 

The time of explosion for ID simulations is quite pro- 
nounced and occurs near the time when the ratio Tgdv/Tq 
becomes of order unity or larger for the last time. How- 
ever, for 2D simulations, all measures show a gradual 
trend toward explosion over 50 to 100 ms timcscales. 
Therefore, to define ioxp we identify the last time that 
Tadv is less than 50 ms (65 ms) for IS-M© (11.2-Mq) 
models. This definition establishes the time of explosion 
after there is a clear trend toward explosion, but before 
post-shock material is launched outward. Acknowledging 
the ambiguity of such a definition, we assign an uncer- 
tainty to tcxp of ±50 ms. 

Having defined toxp, we now define the conditions at 
explosion. For ID simulations, Mexp is defined at 250 

km and at texp- For 2D simulations, Mcxp is defined just 
exterior to the shock (at roughly 300 to 250 km) at the 
time of explosion. Corresponding to the uncertainty in 
iexp, we determine the uncertainty in Mgxp. Figs. iMlfTSl 
and 1161 show a general upward trend of Q, Q/L^^p^, and 
Mgnin around the time of explosion. The values in Table 
[2] are obtained by averaging Q, Q/L^p^, and A/gain for a 
50 ms range centered on tcxp- We also experimented with 
a range of 10 and 100 ms and obtained similar results. 
For ID simulations, we note the conditions at the time 
of explosion, but due to large nonlinear oscillations the 
specific values are sensitive to the definition of toxp- 

FigurefT4lcompares the time evolutions of the net heat- 
ing rates, Q, for the 15.0-lD and 15.0-2D3 sequences. 
While the top panel shows Q for the sequence of ID sim- 
ulations near explosion, the bottom panel shows Q for 
the sequence of 2D simulations (solid) near explosion. 
For comparison, ID models (dotted) for the same range 
of as the 2D sequence are included. The dots show 
Q at the time of explosion for the 2D simulations. As ex- 
pected, Q correlates with L^^ (Fig. fT4|) . Consequently, 
the range of Q for the ID exploding models (~1 to ~11 
B/s, where 1 Bethe (B) = 10^^ erg) is higher than the 
range of Q for the 2D exploding models (^^0.5 to ^6 B/s). 
The ID sequence in the top panel shows pronounced os- 
cillations in Q that correspond to the oscillations in shock 
radius in Fig. [2l The maxima in Q coincide with the min- 
ima in the shock ra dius, which has bee n observed before 
for ID simulations (jBuras "eraI]l2006bD . 

The 2D sequence, on the other hand, shows no such 
oscillations in Q. Rather, the bottom panel of Fig. [14] 
shows a relatively smooth evolution with time. After 



bounce, Q evolves secularly downward for non-exploding 
ID and 2D simulations. In contrast, 50 ms to 100 ms 
before ioxp, Q for exploding 2D simulations remains fiat 
with time or even begins to rise slightly. In either case, 
there is a clear bifurcation in the evolution for exploding 
2D simulations and non-exploding ID and 2D simula- 
tions. 

Similar to Fig. [TH in Fig. [TSjwe plot the heating ef- 
ficiency, Q/Lp^p^. Even after dividing by L^^p^ the effi- 
ciency has a positive correlation with i^^. This suggests 
that the correlation between Q and L^^ is steeper than 
linear. The efficiencies at texp for the ID sequence range 
from 3% to 10%, while the efficiencies of the 2D simula- 
tions at icxp are 3.4%, 3.5%, and 4.4% for L^^ = 1.8, 1.9, 
and 2.0, respectively. Naturally, the trends in Q/L^^p^ 
are similar to the those of Q, including the bifurcation of 
QjLy^p^ for exploding and non-exploding models. 

Fig. [TBI shows the time evolution of Mgain for the 15.0- 
ID and 15.0-2D3 sequences. Once again, the oscillations 
in Afgain correspond to the oscillations in shock radius, 
but this time the minima in Aigain correspond to the 
minima in shock radius. In the exploding models of 15.0- 
2D3, there is a clear secular increase in the amount of 
mass in the gain region leading up to and past explosion. 
This is consistent with the idea that the SASI helps to 
increase the amount of mass in the gain region, tipping 
the scales toward explosive solutions. However, we note 
in Fig. [3|that even the non-exploding models show some 
€ = 1 oscillations, and hence, the SASI is necessary but 
not sufficient in increasing the mass in the gain region. At 
icxp, A/gain = 0.099, 0.0113, and 0.0152 Mq for L^^ = 1.8, 
1.9, and 2.0, respectively. 

We summarize the conditions at explosion in Table 
[H We list for each sequence in Table [Tj the electron- 
neutrino luminosity, L^,^ (column 1), time of explosion, 
toxp (column 2), accretion rate at the time of explosion 
near the shock, A^oxp (column 3), the net heating rate, 
Q (column 4), the heating efficiency, QjLv^u^ (column 
5), and the mass in the gain region at the time of ex- 
plosion, Afgain (column 6). The mass accretion rate at 
icxp is addressed in the next section. The 11.2-lD se- 
quence manifests a wide range of Q, from 0.63 to 4.46 
B/s, and the range of the corresponding heating efficien- 
cies is similarly wide, from 2.4% to 14%. In contrast, 
the 15.0-lD models have much larger Q at toxp, but with 
a smaller range. Typically they are ~3.5 B/s with one 
at 5.57 B/s. The heating efficiencies are similarly high, 
~5.9 to ~10%. On the other hand, 2D simulations have 
lower efficiencies at explosion. All 2D models that use 
the 15-M0 progenitor have Q that lies in the range ^^1 
to ~2 B/s, and have efficiencies that range from ~3 to 
~4%. For most models, Afgain at texp is ^0.01 M0, but 
can be as low as 0.0073 M0 and as high as 0.0319 M©. 

In Figs. [HI [151 and [iBl the exploding 2D models 
show distinctly different values of these indicators com- 
pared to those of the non-exploding ID and 2D simula- 
tions. Q, Q/ij/^j/^, and Aigain trend downward for the 
non-exploding models. Just before, during, and after 
explosion, the heating rate, heating efficiency, and gain 
mass increas e with time for the 2D exploding models. 
IScheck et al.l (|20081 have shown a similar trend in com- 
paring ID and 2D models at the same neutrino luminos- 
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ity, and wc show that this bifurcation exists among the 
2D exploding and non-exploding models as well. 

8. CRITICAL LUMINOSITY AND MASS ACCRETION RATE 

Figure [17] indicates that the concept of a critical lu- 
minosity versus accretion rate condition is relevant for 
time-dependent ID and 2D simulations and quantifies 
the difference in critical luminosity between ID and 2D 
simulations. It plots electron- plus anti-electron-neutrino 
luminosity, ijy^p^ , vs. accretion rate, M . The trajectory 
of each model evolves from right to left as the accre- 
tion rate decreases with time and the luminosity is held 
constant. For models that explode, we plot the lumi- 
nosities and accretion rates at explosion. These points 
strongly suggest L^^y^-M curves that separate exploding 
from non-exploding models. Indeed, models that do not 
explode remain in the low luminosity and high-accretion- 
rate realm that is delineated by the critical curves. With 
error bars, we show the range of accretion rates encom- 
passed by ^oxp ± 50 ms. In general, ID simulations are 
represented by orange, 2D simulations by green hues, 
and 2D-90° runs by blue hues. The three resolutions are 
represented by stars (1), upside-down triangles (2), and 
squares (3). The fi t for L^.p. as a function of M from 
iBurrows fc Gosh^ (|T99l is" plotted (sohd fine). While 
their fit passes through the 15.0-lD points, it overpre- 
dicts the critical luminosities for the 11.2-lD runs by 
~15%. The critical luminosity for all 2D simulations is 
^--^70% of the critical luminosity for ID calculations. 

It is interesting that the points for the 15.0- 
ID sequence ar e con sistent with the results of 
IBurrows fc Goshvl ()1993D . However, we caution against 
interpreting this coincidence too literally. While we do 
employ very similar heating and cooling terms, it is not a 
priori clear how the many differences between our more 
realistic time-depend ent models should comp are to the 
assumptions made in IBurrows fc Goshvl ()1993f ) to obtain 
steady-state solutions. 

Furthermore, we note that the slope of our res ults dif- 
fer from the slope of the [Burrows fc Goshvl (|1993[ ) fit and 
the trend is disconti nuous between the 1 1.2 an d 15 Mq 
results. The fit from IBurrows fc Goshvl (|1993[ ) assumed 
a specific mass for the PNS, while our simulations have 
masses that naturally arise from the initial Fe-core mass 
and the subsequent accretion rates. Just after the Si/0 
interface accretes onto the PNS, the mass of the PNS, 
A'/pNS: is ~1.3 M0 baryonic for the 11.2-M0 progenitor 
and ~1.45 Mq baryonic for the IS-M© progenitor. After- 
wards, the progenitor-dependent accretion rates continue 
to increase the mass until explosion. Of the models from 
the 15.0-2D3 suite that exploded, the PNS masses at ioxp 
are 1.63, 1.59, and 1.52 Mq for L^^ = 1.8, 1.9, and 2.0. 
In contrast, the masses from the 11.2-2D3 suite are 1.4, 
1.38, 1.35, 1.33, and 1.32 Mq for V = 0.8, 0.9, 1.0, 1.1, 
and 1.2, respectively. Assuming that the explosion mech- 
anism can be described by the transition fro m steady- 
state accre tion to a neutrino-driven wind (see iBurrowsl 
([1981 and IBurrows fc Goshvl ([T993h . and work in prepa- 
ration), we expect the critical luminosity to have a mass 

dependence of L^^ oc ^/p^g- Interestingly, this scaling 
with PNS mass resolves the slope discrepancy and the 
discontinuity between the 15 and 11.2-M0 models in Fig. 



Although the concept of a critical and M con- 
dition was derived using steady-state accretion models 
([Burrows fc Goshvlll993l ). Fig. [T7[ shows that the condi- 
tion is appropriate for time-dependent ID and 2D mod- 
els. We find that the critical luminosity for explosion is a 
function o f M, and that it ro u ghly f ollows the trend sug- 
gested by IBurrows fc Goshvl ([19931 ). The critical lumi- 
nosity for explosions in 2D is ~70% of the luminosity re- 
quired for explosions in ID simulations. While there are 
differences between the 90° and 180° simulations and be- 
tween simulations with different resolutions, they roughly 
show the same reduction in critical luminosity when we 
go to 2D. 

9. CONCLUSIONS 

Since IBurrows fc Goshvl ([1993D pubhshed the criti- 
cal luminosity condition for explosions by the neutrino 
mechanism, the relevance of this condition for ID and 
2D time-dependent simulations has remained unresolved. 
Furthermore, recent simulations have hinted that 2D 
simulations teeter on the verge of explosion, when ID 
simulations completely fail. One wonders if the concept 
of a critical condition in 2D simulations can quantita- 
tively explain the trend toward successful explosions by 
the neutrino mechanism. 

With 95 simulations that parameterize L^^, M, reso- 
lution, and dimensionality, we investigate the criteria for 
neutrino-driven explosions. The results of which indicate 
that: 

• A critical luminosity and accretion rate con- 
dition is ob served in time- dependent ID 
simulations. IBurrows fc Goshvl ([l993i ) used the 
success and failure of obtaining steady-state ac- 
cretion solutions to suggest a critical luminos- 
ity and mass accretion rate condition for explo- 
sions. Others have found a critical neutrino lu- 
minosity conditio n for explosion in time-dependent 
ID sim ulations ([Ohnishi et al.l [20061 : iBuras et al.l 
l2006b( ). We note that a critical luminosity is ap- 
plicable for time-dependent ID simulations and 
that it is dependent up on M, as suggested by 
IBurrows fc Goshvl ([l993f ). 

• Radial oscillations are characteristic of 
the transition between ID accretion (non- 
exploding) and exploding simulations. Our 

results add to a growing body of literature that 
obt ain radial osc i Uations near the crit i cal luminos- 
itv (iMavld Il985l : lOhnishi et aT] l2006t iBuras et all 

• A critical L^^ and M condition distinguishes 
accretion and exploding models of time- 
dependent 2D simulations. 

• Unlike ID simulations, radial oscillations are 
not prominent near the critical luminosity of 
2D simulations. Rather non-radial SASI os- 
cillations characterize 2D simulations, even 
for non-exploding models. 

• We suggest that Tadv/Tg is a useful diagnos- 
tic for explosion, but may not be a rigorous 
condition for explosion. Specifically, we found 



that some ID models strongly satisfied this condi- 
tion several times before explosion, and 2D mod- 
els that exploded exhibited trends toward explosion 
well before this condition was met. 

Several integral quantities such as heating 
rate, heating efficiency, and mass in the gain 
region grow with time in ID and 2D ex- 
ploding models, while they de cline in ID 
and 2 D non-exploding models. IScheck et alJ 
compared these quantities for ID and 2D 
simulations with the same neutrino luminosity and 
found declining evolutions in the ID non-exploding 
model and increasing trends in the 2D exploding 
model. We find similar trends. From these re- 
sults, one might suggest that the SASI leads to 
increases in gain region mass and net heating rate, 
which in turn leads to explosion. However, 2D sim- 
ulations that do not explode despite having pro- 
nounced SASI oscillations show declining curves. 
Furthermore, ID exploding models also show in- 
creases in heating rate and gain region mass at ex- 
plosion. Therefore, we conclude that the rise in 
these quantities near explosion are not a unique 
result of the SASI. Instead, the upward evolution 
is correlated with explosion, whether it be in ID 
or 2D simulations. This, combined with the fact 
that the integral measures change from declining to 
increasing evolution near explosion, suggests that 
the upward trends are a consequence of explosion 
rather than a cause of explosion. 

We suggest an improved measure of the res- 
idence time of matter in the gain region, r^, 
that better elucidates the conditions lead- 
ing to explosions. While Tadv is a convenient 
measure of the timescale that a parcel of mass is 
exposed to net heating, we find Tr a more infor- 
mative measure for diagnosing the conditions and 
criteria leading to and during explosion in 2D sim- 
ulations. For example, by following the trajectories 
of particles that start outside and within the gain 
region, we note that much of the mass that accretes 
through the shock receives very little heating as it 
quickly advccts onto the inner core. On the other 
hand, there is a significant amount of mass within 
the gain region that has long residence times and 
receives a great deal of neutrino heating. It is this 
latter material that experiences longer and longer 
residence times and prolonged heating that likely 
reverses accretion into explosion. 

Resolution does affect the SASI and the 
post-shock flow in 2D simulations, but 
hardly alters the critical luminosity and ac- 
cretion rate condition. It noticeably alters the 
explosion time, but the correspondence between 
texp and resolution is not monotonic, nor is it easy 
to disentangle the trends. However, the 2D critical 
luminosity vs. M curve is little affected by resolu- 
tion. 

The 90° simulations explode earlier than 
the 180° counterpart, but again the criti- 
cal luminosity condition is hardly affected. 



We con trast this with the results of iBuras et al.l 
()2006al ) who found an explosion for a 180° simula- 
tion, while a 90° simulation did not explode. How- 
ever, their 90° simulation is ±45° of the equator 
and suppresses £ = 1 and i = 2 nonradial flow, 
the dominant components of the SASI. We simu- 
lated the region from the pole to the equator, which 
suppresses i = 1 but a l lows £ = 2. In agreement 
with our results. iMarekl (j2007l ) found earlier explo- 
sions for 90° compared to 180° simulations using 
th e same hydr o dynam ics and neutrino transport 
as iBuras et al.l ()2006al) but a 90° grid similar to 
ours. In any case, IBuras et al.l ([20063) speculated 
that their 90° simulation was close to explosion and 
not very different than the 180° simulation, and we 
show that the critical luminosity for simulations us- 
ing 90° and 180° are very similar. 

• The critical luminosity for explosions in 2D 
simulations is ^70% of the luminosity re- 
quired for explosions in ID simulations. Ir- 
respective of resolution or the angular size of the 
simulation, this conclusion holds. 

By employing local heating and cooling algorithms, we 
avoid expensive neutrino transport, but we also sacri- 
fice accuracy. In our simulations, heating and cooling by 
neutrinos are decoupled from one another, which is not 
consistent with neutrino transport calculations. In ad - 
dition, simulations using MGFLD iDcssart ct al.' '20& 
and "ray-by-ray" Boltzmann transport ( Marek fc Janb 
|2007() have neutrino luminosities that decrease by ~50% 
from 100 ms to 300 ms past bounce. Our constant L^^ 
certainly does not refiect this evolution in neutrino lumi- 
nosities. Furthermore, a large f raction of is de rived 
from accretion at early times (jScheck et al.l [20061 ). In 
our simulations, L^^ is completely decoupled from such 
effects, and, hence, we ignore any feedback due to this 
coupling. While our assumptions compromise the accu- 
racy of the results, they enable the parameter study that 
we present, and we suspect that the general conclusions 
and trends regarding the critical luminosity condition are 
st ill relevant for simulati ons including detailed transport. 

iMarek fc Jankal l|2007f ) suggest that using a softer EOS 
in conjunction with general relativistic gravity leads to 
more compact PNSs and more favorable conditions for 
explosion; we include neither effect. The soft EOS that 
accompanies explosion has an incompressibility of nu- 
clear matter, K, of 180 MeV, while the stiffer EOS, 
which does not explode, has K = 263 MeV. While 
this suggests an interesting dependence of explosion on 
the EOS, we note that laborator y experiments indi- 
cate that K = 240 ± 2 MeV ((Shlomo et all l2006l : 
iLattimer &: Prakashll2007[ ). a value that is more consis- 
tent with the non-explosive model. Regardless, with a 
compact PNS, the post-shock flow is positioned deeper 
in the potential, making it harder to explode. All else 
being equ al, this should in c rease the critical luminosity. 
However, iMarek fc Jankal ([2007f ) argue that the more 
compact PNS is hotter, emitting more energetic neu- 
trinos (though they are gravitationally redshiftcd) and 
at a higher luminosity, which might compensate for the 
higher critical luminosity. Because we parameterize the 
luminosity and temperature in this study, these effects 
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arc not addressed here and must await a more consistent 
treatment. 

Of the conclusions Usted above, the most striking is 
that time-dependent ID and 2D simulations show a crit- 
ical neutrino luminosity and accretion rate condition for 
successful explosions. Equally significant is that the crit- 
ical luminosity for 2D simulations is ~70% of the critical 
luminosity for ID simulations. This quantitatively sup- 
ports suspicions that multi-dimensional effects increase 
the likelihood of successful explosions by the neutrino 
mechanism. While none have explicitly explored these 
results, the simulations of the past decade are consis- 
tent with this conclusion, while also being intuitive and 
concise. In this sense, these results are expected and 
reassuring. Remarkably, this global condition for suc- 
cessful explosions, which was informed by ID steady- 
state solutions, survives despite the complexities of time- 
dependent multi-dimensional simulations. 

For all permutations of the core-collapse mechanism 
considered in this paper, oscillations are ubiquitous dur- 
ing the transition from accretion to explosion. Interest- 
ingly, the character of these oscillations for ID and 2D 
simulations is quite different, and leads to distinctly dif- 
ferent critical luminosities. Specifically, we suggest that 
the extra degree of freedom in 2D simulations offers an 
easier route to explosion via the SASI. This suggests that 
for other than the lightest progenitors, the core-collapse 
mechanism is inherently multi-dimensional, and an ac- 
curate theory for the core-collapse mechanism must be 
developed in the context of the multi-dimensional Uni- 
verse. 

Our results show that the critical luminosity for 2D 
simulations is lower than the critical luminosity in ID 
simulations, but why is this the case? The most ob- 
vious cause is that the residence time of matter in the 
gain region is different for ID and 2D simulations, and 
the longer residence times in 2D simulations reduce the 
neutrino luminosity required for explosion. So. what ac- 
counts for the different residence times when the accre- 
tion rate is the same? After all, we have shown that 
during the steady-state accretion phase, the advection 
timescales determined by the flow across the gain region 
(T'adv, as calculated by eq. and by = M^^^^/M 

are similar. Put another way, the solid-angle, average, 
radial velocity at radius r is 



Vr{r) 



M{r) 



(15) 



and if p{r) and M{r) are the same for ID and 2D sim- 
ulations, the flow across the gain region should be simi- 
lar. Yet, 2D simulations explode when ID simulations do 
not. This discrepancy can be resolved by noting that the 
flow in the gain region is by no means spherically sym- 
metric in 2D simulations. The plots, histograms, and 
entropy maps of Figs. [8lfT2l show that much of the ma- 
terial that accretes through the shock is channeled to 
plumes and quickly accrete onto the PNS. At the same 
time, there are large regions of high entropy material that 
have long residence times in the gain region. When the 
residence times within these high entropy regions reach 



some threshold, it seems as though sufficient heating has 
occurred to launch an explosion, and because relevant 
regions behind the shock are in sonic contact, the shock 
radius expands globally. Hence, the critical luminosities 
for 2D simulations are lower because some matter resides 
for longer times in the gain region for the same M . 

Since the concept of a critical luminosity condition ap- 
plies for ID and 2D simulations, we assume that it should 
also apply for 3D simulations. Increasing the dimension- 
ality from ID to 2D. has proven favorable for successful 
explosions by the neutrino mechanism. Specifically, the 
extra degree of freedom in 2D simulations increases the 
dwell time of matter in the gain region. Similarly, we 
suspect that differences between 2D and 3D simulations 
will continue the trend toward lower critical luminosities 
for higher dimensions. For example, convection in 2D 
simulations experiences a turbulent energy cascade from 
smaller to larger scales. In Nature and in 3D simulations, 
however, the turbulent energy cascade is described well 
by Kolmogorov theory, in which energy cascades from 
larger to smaller scales. If we analyze the transport of 
matter through the convective gain region as a diffusive 
problem, then the time to diffuse through the gain region 
would roughly scale as 



tdiff 



(16) 



where d is the number of degrees of freedom, I is the 
characteristic size, is the mean-squared speed of the 
diffusing particle, and Tcoh is the average time between 
collisions. The extra degree of freedom in 3D simulations 
alone could increase the dwell time by ^--^50% compared to 
2D simulations. In addition, we suspect that the smaller 
scales of 3D simulations would produce shorter Tcou com- 
pared to 2D simulations and, consequently, further in- 
crease the dwell times. Whether these scaling arguments 
bear out in realistic 3D simulations, in which advection 
and the SASI are key in the flow, will be determined by 
future 3D simulations. If they are relevant, then the in- 
creased dwell time might result in a critical luminosity 
in 3D simulations that is even lower than the critical lu- 
minosity in 2D simulations. In turn, this might lead to 
successful explosions in 3D radiation-hydrodynamic sim- 
ulations. 
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TABLE 1 

Model parameters.-'- 



Sequence Name 


Mass (Mq)'^ 


DimcnsioiT^ 


A^r -R] 


nax (km)° 


ii/e (lO^'^ erg s '-)° 


15.0-lD 


15.0 


ID 


700 


4000 


1.6-2.9 


1 ^ n 9r)1 


15.0 


2D 


250 


4000 


1.5-2.0 


15.0-2D2 


15.0 


2D 


400 


4000 


1.5-2.0 


15.0-2D3 


15.0 


2D 


400 


1000 


1.5-2.0 


15.0-Ql 


15.0 


2D-90° 


250 


4000 


1.5-2.0 


15.0-Q2 


15.0 


2D-90° 


400 


4000 


1.5-2.0 


15.0-Q3 


15.0 


2D-90° 


400 


1000 


1.5-2.0 


11.2-lD 


11.2 


ID 


700 


4000 


1.0-1.8 


11.2-2D1 


11.2 


2D 


250 


4000 


0.7-1.2 


11.2-2D2 


11.2 


2D 


400 


4000 


0.7-1.2 


11.2-2D3 


11.2 


2D 


400 


1000 


0.7-1.2 


11.2-Ql 


11.2 


2D-90° 


250 


4000 


0.7-1.2 


11.2-Q2 


11.2 


2D-90° 


400 


4000 


0.7-1.2 


11.2-Q3 


11.2 


2D-90° 


400 


1000 


0.7-1.2 



This table summarizes the 95 simulations presented in this paper. These simulations represent a 
four-dimensional parameterization that investigates the dependence of the conditions for explosion 
on the accretion rate (column 2), dimensionality (column 3), resolution (columns 4 & 5), and 
neutrino luminosity (column 6).^ Progenitor model."^ Dimensionality.^ Number of radial or effective 
radial zones. ^ Radius of the outer boundary^ The range of neutrino luminosities investigated. 
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TABLE 2 

Conditions at the time of explosion. 



L,, (1052 erg S-l)l iexp (ms)2 Mexp (Mq/s)3 Q (B/s)4 QIU^^J^ Mg^in (Mq)^ 



11.2-lD 



7 



1.3 
1.4 
1.5 
1.6 


932 
496 
312 

409 


0.084_g (,(,3 

112+°°12 

^•11^-0.010 

l36+° °25 
U-iJO_o.oil 

Q 124+0.002 


0.63 
1.50 
2.00 
4.46 


0.0244 
0.0535 
0.0668 
0.1393 


0.0123 
0.0121 
0.0153 
0.0228 


15.0-lD 


2.6 


718 


n 227+0.017 
-0.015 


3.69 


0.0710 


0.0105 


2.7 


459 


U.ZOO_g Q13 


3.79 


0.0701 


0.0129 


2.8 


335 


n O12+0 001 


5.57 


0.0995 


0.0177 


2.9 


216 


n 341+0.027 
U.J4i_Q 027 


3.42 


0.0589 


0.0319 


15.0-2D1 


1.9 


581 


247+0.001 
U.^4'_0.ooi 


1.49 


0.0392 


0.0123 


2.0 


325 


310+0 002 
U.01U_o 001 


1.97 


0.0494 


0.0167 


15.0-2D2 


1.9 


359 


,-, qi O+O.OOl 

U.010_o 023 


1.51 


0.0399 


0.0120 


2.0 


364 


n ni Q+o.ool 

u.olo_o 027 


1.60 


0.0401 


0.0133 


15.0-2D3 


1.8 


778 


2in+0'0i'i 

U.ZiU_o.007 


1.21 


0.0337 


0.0099 


1.9 


649 


n 249+0 001 
u.z4a_Q 015 


1.33 


0.0351 


0.0113 


2.0 


395 


n 299+0 014 
u.zya_o 031 


1.76 


0.0440 


0.0152 


15.0-Ql 


1.7 


837 


202+0'007 
U.ZUZ_o 005 


0.91 


0.0268 


0.0074 


1.8 


419 


n 200+0.029 
u.zoo_o 024 


1.19 


0.0330 


0.0098 


1.9 


418 


n 2S4+0'028 
O.^S4_o.o24 


1.31 


0.0344 


0.0109 


2.0 


361 


n ni Q+o.ool 

U.aici_o 025 


1.71 


0.0428 


0.0152 



15.0-Q2 



1.7 


748 


n 217+0.018 
^'■^'^ ' -0.009 


0.91 


0.0267 


0.0073 


1.8 


517 


n 249+0 011 
U.Z4y_o 002 


1.30 


0.0362 


0.0119 


1.9 


471 


n 2c;s+0 023 

U.ZOO_Q 010 


1.19 


0.0314 


0.0100 


2.0 


299 


„ oifi+O.OOS 
U.JiU_o 001 


1.54 


0.0385 


0.0126 



15.0-Q3 



1.7 


695 


235+0'0l'l 
U.zoo_o 018 


1.14 


0.0335 


0.0096 


1.8 


602 


n 249+0 001 

u.z4y_0 001 


1.08 


0.0301 


0.0089 


1.9 


420 


n 2S9+0 029 
u.zoz_Q 024 


1.42 


0.0373 


0.0122 


2.0 


383 


„ or,7+0.004 
U.OU/ _Q 033 


1.66 


0.0415 


0.0135 



11.2-2D1 



0.7 


808 


n oQi+o ooo 
u.uyi_o 001 


0.23 


0.0165 


0.0050 


0.8 


585 


n nQS+0'005 
0.098_o.oo2 


0.30 


0.0190 


0.0055 


0.9 


471 


n 11C+0.006 

U-118_o 013 


0.43 


0.0242 


0.0068 


1.0 


335 


129+0 020 
u.izy_Q 005 


0.55 


0.0276 


0.0086 


1.1 


288 


n 147+0.016 

^•1^'-0.019 


0.89 


0.0404 


0.0131 


1.2 


288 


n 147+0.016 

^•1^'-0.019 


0.89 


0.0370 


0.0131 



11.2-2D2 



0.7 828 0.092^0 002 0.24 0.0170 0.0054 

0.8 595 0.09710 002 0.31 0.0193 0.0053 

0.9 529 0-104looJ5 0.40 0.0220 0.0070 

1.0 339 0.128lO;oj4 0.54 0.0268 0.0077 

1.1 235 0.163lo;o54 0.85 0.0386 0.0115 

1.2 235 0.163lo;o54 0.85 0.0354 0.0115 



15 



TABLE 2 — Continued 



L.e (10^2 erg s-i)i toxp (ms)2 Mexp {Mq/s)3 Q (B/s)* Q/L,,^,^ Afg,;,, {Mq)^ 



11.2-2D3 



0.8 
0.9 
1.0 


791 
553 
372 


U.Uyi_Q QQJ^ 

101+0009 
n intr + O.OOS 


0.32 
0.40 
0.60 


0.0200 
0.0225 
0.0300 


0.0057 
0.0059 
0.0087 


1.1 


257 


U.iDZ_Q Q25 


0.91 


0.0413 


0.0122 


1.2 


257 




0.91 


0.0378 


0.0122 



11.2-Ql 



0.7 


777 


091+0 001 


0.23 


0.0165 


0.0048 


0.8 


500 


Q 11 1+0.012 


0.31 


0.0191 


0.0060 


0.9 


784 


091+0 001 


0.78 


0.0434 


0.0133 


1.0 


294 


144+0 019 


0.52 


0.0260 


0.0076 


1.1 


177 


U.iDD_Q QQ4 


0.91 


0.0414 


0.0124 


1.2 


177 


166+°°3* 

U.iDD_Q 334 


0.91 


0.0379 


0.0124 



11.2-Q2 



0.7 


631 


ogfi+^ ooa 

U.UyD_(j QQ2 


0.22 


0.0158 


0.0040 


0.8 


525 


n int;+o.oi3 

U.1U0_Q QQg 


0.28 


0.0174 


0.0056 


0.9 


413 


124+0.001 
U-1Z4_0.004 


0.39 


0.0218 


0.0071 


1.0 


293 


145+0018 


0.54 


0.0269 


0.0075 


1.1 


270 


0.158l°;~^ 


0.82 


0.0371 


0.0107 


1.2 


270 


0.158l»:00^ 


0.82 


0.0340 


0.0107 



11.2-Q3 



0.7 


759 


0.09l±»;'»i 


0.22 


0.0160 


0.0047 


0.8 


583 


0.098iHS^ 


0.30 


0.0186 


0.0056 


0.9 


417 


„ 124+0.001 


0.39 


0.0218 


0.0065 


1.0 


293 


145+0018 


0.56 


0.0282 


0.0083 


1.1 


256 




0.84 


0.0383 


0.0125 


1.2 


256 


U.iDZ_Q Q24 


0.84 


0.0351 


0.0125 



Electron-neutrino luminosity The time of explosion The mass accretion rate at explosion, The net heating rate The heating efficiency, 
where Lj^^j^-^ is the combined electron- and anti-electron-neutrino luminosity.^ The mass in the gain region.^ The table is divided into sections 
corresponding to the conditions at explosion for each sequence presented in Table fTI 
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0.0 0.2 0.4 0.6 0.8 1.0 
Time after bounce [s] 

Fig. 1. — Accretion rate, At, vs. postbounce time above the stalled shock (2.50 km). Th e solid and dashed l ines show the time-dependent 
accretion rate for the 11.2- and 15.0-Mq models, respectively l lWooslev et al . 2002; Wooslev fc Weaveilll995l ). While the outer portions of 
the Fe core accrete {t = 0-50 ms), M is as high as 10 Mq/s and decreases to 2 Mq/s. After the Fe core fully accretes, the accretion rates 
for the two models diverge. For the 11.2-Mq (I5-M0) model, it takes 50 ms (100 ms) to accrete the Si-burning shell. As the Si/O interface 
accretes, M plummets to 0.3 Mq/s (0.2 Mq/s) for the 15 Mq (11.2 Mq) model. Afterward, M declines slowly to 0.2 Mq/s (0.08 Mq/s). 
Together, these two models slowly sweep through a range of accretion rates from 0.3 Mq/s to 0.08 Mq/s, which enables a parameterization 
in M. 
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Fig. 2. — Shock radius, -Rshocki vs. postbounce time for the ID sequences. The top panel displays the radii for 15.0-lD, and the bottom 
panel shows the radii for 11.2-lD. Each line is labeled by the electron-neutrino luminosity in units of 10^^ ergs s~^. Comparing models 
within a sequence, all models show similar behavior prior to the accretion of the Si/O interface. The accretion of the Si/O interface either 
initiates an explosion or excites radial oscillations in the shock radius. For the 15.0-lD sequence, the oscillation periods range from ~90 ms 
for Li,^ = 2.2 to ~170 ms for L^^ = 2.8. Of the models that oscillate, the lower luminosity simulations decay in the oscillation amplitude, 
while higher luminosity models oscillate and explode. The timescales for decay range from ~450 ms for Li^^ = 2.2 to 1.0 s for L^^ = 2.5. 
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Fig. 3. — Shock morphology vs. postbounce time for the 15.0-2D3 sequence. The top panel shows oq or (-f?shock)i the average shock 
radius, and the middle and bottom panels plot the 1=1 and 1 = 2 components divided by to the average shock radius; a^/aQ and 02/00 1 
respectively. The filled circles in the top panel indicate texp, the time of explosion, and the error bars show our estimate of uncertainty 
(±50 ms). Unlike the ID models (Fig. [2}, the 2D simulations show very little radial oscillation near the critical luminosity, but significant 
nonradial oscillations. In general, ai/ao grows from 0% to ~10% for all models, and 02/00 grows to ~5%. The oi/oq amplitudes hint at 
a correlation with Ly^ (see [|4]for a discussion). Near explosion the amplitude of oi/oq can peak around ^^20%. 
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Fig. 4. — The shock morphology vs. time for three different resolutions of 2D simulations (15.0-2D1, 15.0-2D2, and 15.0-2D3). The 
layout is similar to Fig. \3\ The electron-neutrino luminosity for each is 1.9 X 10^^ ergs s^^. From lowest resolution (15.0-2D1) to the 
highest resolution (15.0-2D3), the explosion times arc 581, 359, and 649 ms. Hence, t^xp is not monotonic with resolution. See [Ofor a 
discussion. 
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Fig. 5. — (^Jshock) a-nd O'i/o.o vs. time for two models of 15.0-Q3 and 15.0-2D3. This compares the differences among the 2D runs using 
180° (soUd) and 90° (dashed). L^^ = 1.8 corresponds to the purple curves and L^^ = 1.9 corresponds to the green curves. Inspection of 
(^^shock) (top panel of Fig. [Sjl shows that the 90° simulations consistently explode before the 180° simulations. Other than the time near 
~300 to ~400 ms, all models show similar evolution in the amplitude of 02 /ao. The one exception appears around ~300 to ^^400 ms, when 
the 90°, Li/^ =1.9 model explodes and shows the largest amplitude. 
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Fig. 6. — The ratio of advection and heating timescales, Tadv/^q, vs. postbounce time for the 15.0-lD sequence. The top panel shows 
the models that do not explode by 1.3 s after bounce, and the bottom panel shows those that do explode. Despite dramatic oscillations, 
Tadv/'^q rarely reaches a value of above ~1 for the non-exploding models (top panel). In all non-exploding cases, except L,^^ = 2.5 (in units 
of 10®^ erg s"'^), the oscillations decay (see text and the caption of Fig. [2]for the timescales). For exploding models (bottom panel), r^^^/rq 
makes large excursions from ~0.1 to ~100. Larger luminosities produce more excursions, with one at L^^ = 2.9 and four at Ly^ = 2.9. We 
conclude that T^d-^/Tq > 1 is a useful diagnostic, but not a rigorous condition for explosion. 
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Fig. 7. — Similar to Fig. |6] but we plot r^^^jrq vs. postbouncc time for the 15.0-2D3 sequence. The models that don't explode, 
Lu^ = 1.5, 1.6, and 1.7 (in units of lO^'^ erg s~^) show relatively constant ratios for all times. On the other hand, the exploding models, 
L^^ = 1.8, 1.9, and 2.0, start at low ratios (~0.1) and increase to ratios of order ~10. Unlike the ID models, once this ratio reaches 
large values, it remains there for the rest of the simulation. The times of explosion are 395, 648, and 778 ms for Lj^^ = 1.8, 1.9, and 2.0, 
respectively. Notice that the ratio of timescales begin to rise before T^^^/rq = 1, suggesting that these models are on their way to explosion 
before the condition r^^^/rq > 1 is met. 
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Fig. 8. — The mean residence time {{rr)) (top panel) and the fraction (fr) of particles with Tr > 40 ms (bottom panel) as a function of 
time for iistart = 400 km (solid line) and iJstart = 150 km (dashed line). These residence times were calculated for the 15.0-2D3, = 1.9 
model. 50,000 tracer particles are initiated randomly in solid angle at either -Rstart = 150 km or -Rstart = 400 km. The trajectories for 
these particles are integrated for 150 ms, and a new generation is initiated every 10 ms. The time on the horizontal axis corresponds to 
the start time for a generation. From ~200 ms and beyond, particles that start in the gain region (Hstart = 150 km) consistently show 
larger (rr) and more particles with longer residence times. There are four phases in this plot that correspond to very short residence ti mes, 
steady grow th in the SASI, an increase in the vigor of the SASI, and the initial stages of explosion. See [|6]for a discussion and Figs. 1101 
111! and 1121 for entropy maps. 
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Fig. 9. — The distribution of residence times, Tr, for four generations of tracer particles corresponding to the four phases in Fig. [8] 
The times shown in the legend correspond to the starting time for each generation, and iJstart indicates the starting radius for all 50,000 
particles of each generation. Rstart = 400 km is situated well outside the shock radius, and -Rstart = 150 km is approximately the middle 
of the gain region. The top panel shows that most of the mass that accretes through the shock traverses the gain region on relatively short 
timescales, ~40 ms, while the bottom panel shows that in addition to rapidly accreting plumes, there are regions in the gain region in 
which mass has a large residence time and is subjected to prolonged heating. 
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Fig. 10. — Colormaps of entropy (/cB/baryon) for 15.0-2D3, Ly^ = 1.9 and for the four phases shown in Fig. [9] Each panel is 400 km X 
800 km. These snapshots of the flow correspond to the times in Fig. |9l which label each histogram by the start time of the four generations 
of particles shown. The four phases are: 1) there is very little convection or SASI and the particles accrete quickly {~9 ms) onto the PNS 
(top left), 2) the SASI undergoes a steady growth and the particles initiated at 400 km (200 km) have a mean residence time of 22 ms (32 
ms) (top right), 3) the SASI is becoming quite vigorous and the residence times for both sets of particles are steadily increasing (bottom 
left), and 4) The SASI is very vigorous and 150 ms later (see Fig. I12I I the model undergoes explosion (bottom right). 




Fig. 11. — The same as Fig. 1101 except the flow is shown 75 ms after the start of each generation of particles, which is the midpoint of 
particle tracking. Low entropy plumes that quickly transport newly-accreted material to the PNS are quite pronounced in the bottom-right 
panel. 




Fig. 12. — The same as Fig. 1101 except the flow is shown 150 ms after the start of each generation of particles, which is the end of 
particle tracking. Note that the simulation is exploding in the bottom-right panel. 
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Fig. 13. — Entropy vs. radius for the times shown in Fig. 1101 Solid red lines are entropy profiles for the 15.0-lD, Ly^ = 1.9 simulation, 
and points show the entropy and radius for each zone of th e 15 .0-2D3, Li,^ =1.9 model. These entropy profiles are coeval with the residence 
time distributions in Fig. |5]and the entropy maps in Fig. 1101 Features of the entropy profiles highlight the shock (~150 km for t = 0.130 
s and ~200 km for all other times), the gain region, (~100 km to ~ 200 km), and the cooUng region (~50 km to ~100 km). ID profiles 
show a negative entropy gradient in the gain region that results in convection in 2D simulations. When the SASl dominates the post-shock 
flow, t = 0.370, 0.530, and 0.600 s, the shock is distinctly aspherical, and a large range of entropies characterize the gain region. Low 
entropies correspond to plumes that funnel matter onto the PNS, while high entropies correspond to regions of long residence times and 
more integrated net heating. This is the likely cause for 2D simulations having lower critical luminosities than ID simulations. 
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Fig. 14. — Time evolution of the net heating rate, Q, for the 15.0-lD and 15.0-2D3 sequences. The top panel plots Q for the ID sequence 
only for the range of L^,^ near the ID critical luminosity. The bottom panel highlights the range of near the critical luminosity of 
the 2D sequence, and plots both ID and 2D simulations. The dots show the Q at the time of explosion for the 2D simulations. Q shows 
pronounced oscillations for the ID simulations near explosion (top panel), with the peaks corresponding to dips in shock radius in Fig. 
[2] For non-exploding ID and 2D simulations Q evolves downward. Exploding 2D simulations, however, have Q that tends toward fiat or 
upward evolutions. 
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Fig. 15. — Similar to Fig. 1141 but here we plot the heating efficiency, Q/Ln^i,^ , vs. postbouncc time, 
sequence near the critical luminosity range from 3% to 10%, while the efficiencies of the 2D simulations at texp are 3.4%, 3.5%, and 4.4% 
for Li^^ = 1.8, 1.9, and 2.0, respectively. For non-exploding ID and 2D simulations Q/ Liy^a^ evolves downward. Exploding 2D simulations, 
however, have efficiencies that tend to stay flat or increase. 
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Fig. 16. — Similar to Fig. [14] but here we plot the niciss in the gfiin region, ^Igaim vs. postbounce time. The oscilljitions of A^gg^^j^ in 
the ID models (top panel) correspond to the oscillations in shock radius. The minima in A/gain correspond to the minima in shock radius. 
In the exploding models of 15.0-2D3, there is a secular increase in the amount of mass in the gain region leading up to and past explosion 
that is distinct from the non-exploding ID and 2D models. 
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Fig. 17. — The electron- and anti-electron-ncutrino luminosity, L^^j}^ , vs. mass accretion rate, M, of explosion. The accretion rate is 
defined just exterior to the shock at fexp- Because the time of explosion is difficult to define, we also show the range of accretion rates 50 
ms before and after iexp as error bars. In general, ID simulations are represented by orange, 2D simulations by green hues, and 2D-90° 
runs by blue hues. The points that lie in the M range from ~0.8 to ~1.7 Mq/s correspond to simulations using the 11.2-Mq progenitor, 
and the points in the range from ~0.2 to ~3.5 Mq/s correspond to 15 Mq. The three resolutions are represented by stars (1), upside-down 
triangles (2), and squares (3). The fit for L^^v^ as a function of M from lBurrows fc Goshvl II1993I ) is plotted (solid line). While this analytic 
expression passes through the 15.0-lD points it over predicts the critical luminosities for the 11.2-lD runs by ~15%. These results show 
that time-dependent ID and 2D simulations reproduce the critical L^^ and M condition, the critical luminosity is indeed a function of M, 
and the critical luminosity for 2D simulations is ~70% the critical luminosity for ID calculations. 



